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ABSTRACT 


The report addresses itself to a summary of the numerical study of the 
Reynolds-stress turbulence closure for separating, reattaching, recirculating 
and redeveloping flow. The calculations were made for two different closure 
models of pressure - strain correlation. The results were compared with the 
experimental data. Furthermore, these results were compared with the compu- 
tations made by using the one-layer and three-layer treatment of k-e turbu- 
lence model which were developed iri the first part of this project sponsored 
by NASA-Ames. Generally the computations by the Reynolds-stress model show 
better results than those by the k-e model, in particular, some improve- 
ment was noticed in the redeveloping region of the separating and reattachiny 
flow in a pipe with sudden expansion. 
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coefficients in the R.S. model of turbulence 

<0 

coefficients in the turbulence model 

diameter of pipe downstream of expansion 
diameter of pipe upstream of expansion 
step height (= (D-d)/2) 

—2 

turbulent kinetic energy (= u^/2) 

Nusselt number 
production rates In RSM 
turbulence energy generation rate 
pressure 

redistribution terms in RSM 

wall redistribution terms in RSM 

Reynolds number based on diameter of pipe downstream 

of expansion 

radial coordinate 

mean velocity in x direction 

turbulent fluctuating velocity components 

mean velocity in r direction 

coordinate parallel to pipe axis 

nearest distance of the node from the wall 

effective diffusivity 

dissipation rate of turbulence energy 

dynamic viscosity 

effective viscosity (= p + p t ) 


V 
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V 


a 
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Subscripts 

1 ,J 

z 


turbulent dynamic viscosity (= C pk /e) 

kinematic viscosity 

density 

Prandtl number 

turbulent Prandtl numbers for diffusion of k and 
shear stress 
dependent variable 


tensor notation 

turbulent values of quantity 


1. INTRODUCTION 


The heat transfer augmentation which occurs as a result of flow separa- 
tion and subsequent reattachment is a very important aspect when one is com- 
puting the flow over an aerodynamic wing w?th separation bubbles. Wh 1 ’ le there 
have been many significant contributions to the flow of this type, current 
understanding of the process of these flows is still relatively poor partly 
because the flow mechanism is complex and partly because presently existing 
turbulence models still have many limitations for a prediction of a wide range 
of parameters in the separating, reattaching and recirculating flows. 

The effect of the recirculating flow was examined by Tani et al. [Ij re- 
porting that the maximum back flow velocity is usually over 20% of the free- 
stream velocity. Bradshaw and Wong [2] and Smyth [3] studied the character- 
istics of the flow near the reattachment point and found that the turbulence 
characteristics in the shear layer are usually transported to the downstream 
of the reattachment. 

The above experimental observations suggest that an elliptic approach is 
always required for reattaching flows. Moreover, it suggests that the large 
scale-eddies that are developed in the separated free-shear layer persist in 
the reattachment and redeveloping regions. It was also concluded that the 
turbulence models which are applicable for mixing-layer flows are also ap- 
plicable for the reattaching shear flows. In fact, many engineers and re- 
searchers have been computing these flows by adopting the standard k~e 
model because of its simplicity and success in the prediction of the flow 
field. However the predictions by the standard k~e model are not always 
successful for higher Reynolds number flows or for larger step heights. This 
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is because the turbulence shear stresses in the separating flow are not com- 
puted correctly by using the k~ e model. 

Despite the limitations of the k„e model, the prediction 'of the 
turbulence quantities could be improved by employing near-wall models of 
Chi eng and Launder [4] and Amano [5 and 6] in which the local variations of 
turbulence quantities are evaluated through the viscous sublayer and the 
overlap layer regions. 

A recently simplified version of Reynolds stress turbulence closure was 
developed and successfully applied for separating and recirculating flows by 
Hanjalic, Launder and Sindir (shown in Driver and Seegmiller [7]). This mod- 
el, the so called algebraic-stress model, generally gives better results in 
the recirculating region, however, the results obtained in the region down- 
stream of the reattachment by using the algebraic-stress model were not as 
good as the prediction by the standard k~e model. 

In this paper the full Reynolds-stress transport equations are reviewed 
and adapted tor the elliptic flows. For the pressure-strain correlation, two 
different formulations were applied; one is proposed by Naot et al. [8], and 
the other by Launder et al. [9]. Because of the complexity of these equations 
some significant modi fi cat .ons had to be made in the corresponding momentum 
equations in evaluating diffusive terms by employing the Boussinesq viscosity 
concept. These modifications are described in the following section. The 
result of the calculations are presented in the last major section and are 
compared with the experimental data of the authors* previous work [6]. Com- 
parisons are also made to the computational results by the k~ e model 
with different near-wall models. 
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2. MATHEMATICAL FORMULATION 

2.1 Equations and Turbulence Models 

The governing equations for the models employed in this study can be 
written in the following general form 

7 [fj (rpU+) + (rpV*)3 

■ ? & # + If < rr 2 U) 

All the equations solved are summarized in Table 1. 

Besides the conventional form of the standard k~e model, we con- 
sider three different turbulence models. 

(Model 1) Boussinesq Viscosity Model 

This is the ordinary k~e model, but the Reynolds stresses are com- 
puted separately so that this model can be replaced with the Reynolds-stress 
model (RSM) easily. By using the Boussinesq viscosity concept, all the shear 
stresses may be expressed as follows: 


pu z * 2/3 pk-2ji t !£, 


pv 2 * 2/3 pk-2v t |i, 


pw 2 = 2/3 pk-t2pj. V/r}, 


( 2 ) 


pUV = - 
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The RSM consists of transport equations for four stresses (u , v , 

~~2 

w and uv) and each term corresponding to the one in Eq. (1) is listed in 
Table 1. In these equations the productions for these stresses are given as 


P„ « . 2p (u 2 & 

11 M v ax ar' 


2 aV 


aV, 


P 22 = - 2 p (v fc H + uv 


ar 


3X ‘ 


(3) 




P 33 = -{2pw V/r) 


p i2 = " p (y2 If + u2 lx ■ {uv v/r}) * 


In this study two different redistribution terms (pressure-strain corre- 
lation) are used. One is proposed by Naat et al. [ 8 ] and the other by Launder 
et al. [9]. 


(Model 2) KSM with redistribution of [ 8 ] 


R 11 = ~ C «j> 1 pgI< " C 4>2^ P 1 1 ‘ 3 P )* 


R 22 s " C 4 >l pek ^ y2 " " C <j>2^ P 22 ' 3 P )* 


R 33 ~ “ C 4» l pe{< ’ C $2^ P 33 3 P ^ 


(4) 


-1 


‘12 


= - C^pek UV - C^P 12 . 
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(Model 3) RSM with redistribution of [9] 


R n ■ ' c +i» ek " , (“ ? • fo - B i< p n -! p J 


2B 2 »4 + 2B 3 (p7^ + puv|^ + ^ ! 


H Zi ■ - Sl^'V - §0 -B,(P 22 -fp) 


- 2B Z> k J? + 2 h^ Z sr + M fr + 7 P >- 


(5) 


• 1/...2 2 , 


R 33 = " C il pck ( w “ 3^ ” B 1^ P 33 “ 1 P ^ 


- {2B 2 pkV/r} + 2B 3 ( {^} + ]p) 


-1 


3U ^ 3 V ■ 


*12 S - C *l pek uv - B 1 P 12 “ h*®*® 


+ B~p(u 2 ill + v 2 — +'uv [— + — ]). 
3 v sr 3X l 3x 3r ,J/ 


The generation term P used in Models 1, 2 and 3 is given as 


P = - pCuv l|~ + fjf) + fx + y2 I? + iw 2 V/r}). 


(«) 


The constants used in these models are as follows 


Ci C 9 C e C 
M 1 2 S 4 


C <|>2 B 1 B 2 B 3 


0.09 1.44 1.92 0.25 1.50 0.40 0.764 0.182 0.109 
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2.2 Low-Keynolds Number Form of RSM 

For the incorporation of wall effect into RSM, the wall redistribution 
term proposed by Launder et al. [9] is used for the low-ReynoTds number form 
of RSM. The terms of low-Reynolds number form are denoted as R wl ^, R w22 » 

R w33 , and R for u , v , w , and uv respectively (in source term 

in Table 1). These are expressed as follows for the present flow regions. 

«wn ■ c°- 125 f - ! k ) - °- 015 ^ ^ 

. 3/2 

R w 22 ■ [0.125 f (v 2 - 2 k) -0.015 (2uv 

( 7 ) ’ 

R W 33 ‘ °’ 125 F <7 - |k) ^ /2 

<W ’ C0.125 f uv - 0.015 1(7 - 7)lf - |i))]|| 3/2 

It should be noted here that in the above expression of R w ]*|» R W 22» 
and R w i 2 there is a term -0.015 { }. In the original paper, this term 
is positive. The negative sign has been given since the flow geometry here is 

different from that in the original paper and so the velocity gradients had to 

be accordingly converted by changing the sign. 

For the high-Reynolds number form of RSM, these terms, R w ^j» are simply 
set equal to zero. 

2.3 Numerical Method 

Formulation and discretization of Eq. (1) was performed by using the con- 
ventional control volume approach of Patankar [10] by breaking each of the 
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aquations Into diffusive, convective and source terms, The systems of equa- 
tions were made tridiagonai such that they could be solved Iteratively by the 
Tri-Diagonal Matrix Algorithm (T.D.M.A.). 

Laminar viscous diffusive terms had to be incorporated In the main gov- 
erning equations to take Into account laminar viscous terms which are signif- 
icant and adjacent to solid wall boundaries. 

The Reynolds stresses were each programmed into individual subroutines. 
Furthermore, each Reynolds stress has two separate subroutines -- one with the 
simple equation of Model 1, and the other with the transport equation incor- 
porating both Models 2 and 3, 

For the computation of Reynolds stress transport equations, u^Uj is 
contained in diffusive and source terms of its own ujtT. transport equa- 
tions. From the numerical point of view a convergence of these transport 
equations would hardly be attained due to their explicit form. In fact the 
chance of convergence was minimal when the computation started with an ini- 
tially guessed uj'uj field. For this reason we adopted a three-pass pro- 
cedure. Initially the general U and 1/ momentum equations were solved along 
with the conventional k-e equation with the generation. 


" ■ + 


^) 2 + 

ax'' 


+ 

v ax' 


2 ( f ^) 2 + 2 (-“) 2 ] 


( 8 ) 


After about 150 iterations, the simple Reynolds stress equations (Model 1) 
were brought in and the generation term of Eq. (8) was replaced with Eq. (6) 
to incorporate the Reynolds stresses. After about 20 ~ 50 iterations 
therefrom, the simple model (Model 1) was replaced by the transport equations 
of the Reynolds stresses. Complete convergence was achieved after about 
100 - 150 iterations more. 
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For the computation of diffusive terms of RSM, these are divided into two 
groups as 


OU) * i [fj (r P kC s u 2 /c |i) + (rpkc//, ffj j t U'u) 


(*) 


where 


2 2 2 — 
4> 8 u , v , w or uv 


110 ) 


0'(*) » p Cjp (rp. s uv/e |£) + Ip (pkC s uv/c f£)] 


In Eq. (9) the first term is the diffusive term of Eq. (1) whereas the second 
term D'($) is taken in as the source term. This modification was done pri- 
marily to make programming easier and also to make the equation compatible 
with the standard form of the transport equation (Ea, (!}}• 

The ceM structure employed in this study is shown in Fig. 1. The scalar 
quarsUtie* such as p, k, and e are evaluated at grid point P, whereas U and 
V rnomentums are computed at staggered grids. For Reynolds stresses, the shear 
stress uv is evaluated at the left-bottom corner of the scalar cell as is done 
by Pope and Whitelaw [11]. However, normal stresses are computed at different 
points in the manner indicated in Fig. 1. This arrangement has more advan- 
tages in evaluating the terms with shear stress and shear strain combinations 
than those evaluated at the same scalar node point. 

2.4 Boundary Conditions of RSM 

Boundary values of the Reynolds stresses at solid wall have been taken 
from Launder et al . [9], where near the wall 
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UV = - U r 2 + 7 0 

u 2 = 5.1 U 2 
r 

V 2 " = 1.0 U 2 

t 

w 2 a 2.3 u 2 

t 

and since 


original page ® 
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k a 3.5 U 2 , 

T ’ 

the boundary values can be expressed in terms of k as follows 
uv = - 0.286 k + ^ 

u 2 = 1.457 k 
v 2 « 0.286k 


(iu 


(12) 


( 13 ) 


w 2 = 0.657 k 

So along the solid walls, the above values have been prescribed where at the 
exit and the symmetry line, zero gradient is assumed as the boundary condi- 
tions. 

3. PRESENTATION OF RESULTS AND DISCUSSION 

As explained earlier, the computational method employed was a three-pass 
procedure whereas in the conventional k-e model it was solved up to partial 
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convergence and then the Reynolds stresses equations were subsequently brought 
in. Complete convergence for a 32 x 22 grid took about 350 iterations from 
the start with a total central processor time of about 20 minutes on a (JNI VAC 
1100 computer. The grid independency tests were also performed for several 
different non-uniform grid systems and the conclusion was that the present 
system of 1.03 x 0.95 would provide a compromise between numerical stability 
and accuracy. 

The numerical results of the present computational method are compared to 
the experimental data of Amano et al . [6]. The Nusselt number or the con- 
vective heat transfer coefficient was evaluated by employing the interlinkage 
between the near-wall variation and the local wall fluxes developed by 
Jayatilleke and used in Ref. [5]. 

The results presented and discussed here are for a diameter ratio (d/D) 
of 0.586 and Reynold numbers (Re^) of 10,000 and 20,000. Nusselt number 
distribution along the pipe wall downstream from the expansion as a function 
of dimensionless axial pipe length is presented and discussed for different 
models and options of turbulence. Comparison between the computed results of 
these models and experimental data is also made. 

Figure 2 shows the Nusselt number distribution for Reynolds numbers of 
10,000, and 20,000. Computational results are obtained using a high Reynolds 
number form of RSM for both Model 2 and Model 3. The Nusselt numbers computed 
by Model 3 show generally 4 to 6% lower than those by Model 2. In Fig. 2(a) 
(Re n = 10,000), the discrepancy in peak Nusselt number (Nu ) between 
experimental results and those of Model 3 is about -14% while it is about 
-8.5% for Model 2. A similar trend is observed in Fig. 2{b) (Re^ = 20,000), 
the discrepancy in Nu [J|ax between experimental results and those of Model 3 
being -24% and -22% for Model 2. The Nusselt numbers are, however, more 
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elevated an account of higher heat trnsfer effects due to increased velocity 
for Kdp = 20,000. 

The general trend compared to the experimental results is similar, Tne 
Nusselt number starts out at a low value, begins to increase, reaches a peak 
value somewhere 5-8 step heights (x/H) downstream of the expansion ana 
then begins to drop down again to reach a fully developed value. 

Although it is difficult to conclude which model predicts Nusselt number 
better, the results by using Model 3 seem to be lower than those by Model 2 in 
the redeveloping region which demonstrates that the results by Model 3 are 
slightly closer to the experimental data. 

Figure 3 presents the Nusselt number distribution for Rep = 10,000 and 

20.000 with the low-Reynolds number form by using both Model 2 and Model 3. 

For Re h a 10,000, the discrepancy in Nu mav for Model 3 is - 9 % and it is 
- 3 % for Model 2. A similar trend can be seen for Rep = 20,000. The dis- 
crepancy in Nu x for Model 3 is - 9 % and this is - 5 % for Model 2 (Fig. 3(b)). 

In general, although Model 2 predicts a better agreement insofar as 
Nu max is concerned the overall agreement of Model 3 with the experimental 
results is better. 

Figure 4 shows the Nusselt number distribution for Reynolds number of 

10.000 and 20,000 obtained for Model 3 by using both lowand high-Reynolds 
number forms. It is shown that the levels of Nusselt number are higher by 
using low-Reynolds number form than the results by high-Reynolds number form. 
The agreement between computational results and experimental data are general- 
ly good in the reattachment region, however, the prediction by low-Reynolds 
number in the redeveloping region show slightly too high values than that by 
high-Reynolds number model. It has been demonstrated by Launder et al. [9 J 
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that their low-Reynolds number model predicts well for fully developed channel 
flows, However, shear strains are not large enough in the reattaching and 
redeveloping regions. In order to be compatible with the fully developed flow 
the shear stress - shear strain term (the second term in Eq. ( 7 ) ) should be 
sensitized more for the reattaching and redeveloping flows. 

Finally, comparisons between the RSM and the k~e model are shown in 
Fig. 5 for both Re^ = 10,000 and 20,000. For the computation of k-c 
model both the one-layer model and the three-layer model were employed. The 
one-layer model is the near-wall treatment in which the local variations in 
viscous effect are not considered at all, whereas the three-layer model is the 
treatment in which the local variations in viscous sublayer, buffer layer, and 
fully turbulent region are taken into account [5 and 6]. Figure 5 shows that 
the RSM gives the best prediction of all, both in reattaching and redeveloping 
regions. 

4. CONCLUSION 

From this study the following conclusions can be drawn. 

1. The computation of the reattaching and redeveloping flows 
was improved by incorporating the RMS. 

2. The results by Model 2 and Model 3 show about 4-6% differ- 
ence and particularly in the redeveloping region Model 3 
seems to give better results than Model 2. 

3. In the reattaching and redeveloping regions the prediction 
by the low-Reynolds number cannot be improved. It is sug- 
gested that the redistribution term of wall effect needs 
revising for the reattaching and redeveloping flows. 
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4. Incorporation of RSM improves the prediction more than im- 
provement of k*.e near-wall model. It suggests that 
one should employ the RSM model rather than refining the 
k~e model . 

Finally, it could be concluded that, although the RSM improves the pre- 
diction in the separating and reattaching flows, the RSM in its present model 
needs to be modified considerably for the computation of these flows. 
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Table 1 Summary of Equations Solved 
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FIGURE CAPTIONS 


FIG. 1 The grid, storage locations and control volumes. 

FIG. 2 Nusselt number distribution along the larger pipe wall 

a) Re D = 10,000 

b) Re p = 20,000 

FIG. 3 Nusselt number distribution along the larger pipe wall 

a) Re D = 10,000 

b) Re D = 20,000 

FIG. 4 Nusselt number distribution along the larger pipe wall 

a) Re p = 10,000 

b) Re g = 20,000 

FIG. 5 Nusselt number distribution along the larger pipe wall 

a ) Re g - 1 0,000 

b) Re p = 20,000 
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Fig. 1 The grid, storage locations and 
control volumes. • , grid point, 
, control surface. 
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Fig. 2b 
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EXPERIMENT (Amano et al.) 
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